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Abstract 


A  model  for  numerical  simulation  of  the  solar  pumping  of  a  simple  atomic 
system  is  constructed  and  compered  with  the  exact  analytical  solution.  The 
process  is  then  extended  to  cover  the  more  complex  5-level  5 -transition  barium 
ion  term  scheme,  and  the  61-level  86-transition  barium  neutral  term  scheme. 
An  advantage  of  the  stepwise  simulation  is  that  in  addition  to  yielding  the  equi¬ 
librium  relative  level  populations  and  transition  intensities,  it  also  permits  the 
dynamic  grow -in  to  equilibrium  to  be  studied.  The  neutral  barium  system  has 
also  been  studied  with  the  inclusion  of  photoionization  from  each  of  several 
metastable  levels. 
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Photo-Equilibrium  of  Barium 


1.  INTRODUCTION 


When  a  cloud  of  atom  (or  ion)  vapor  is  irradiated  by  sunlight,  the  relative 
populations  of  the  energy  levels  and  the  relative  intensities  of  the  optical  transi¬ 
tions  are  determined  by:  (1)  the  oscillator  strengths  of  the  transitions,  and 
(2)  the  relative  intensities  of  sunlight  at  the  wavelengths  of  these  transitions. 

In  the  case  of  barium  ion,  there  are  5  levels  and  5  permitted  transitions;  the 
equilibrium  solution  when  all  transitions  have  negligible  optical  thickness,  has 
already  been  given  by  Best  (1968). 

If  we  number  the  energy  levels  of  a  term  scheme  from  1  to  W,  starting  with 
the  ground  state,  then  the  rate  of  change  of  population  Pn  of  leve.  i  due  to  absorp¬ 
tion  of  solar  photons  by  lower-lying  energy  levels  is  equal  to 


dP 

n 

dt 


(1) 


where  P  is  the  population  of  the  m1*1  level,  0  is  the  solar  irradiance  (photons 

•9  “  1  1  uiTi 

cm  6  sec"1  Hz"1)  at  the  wavelength  k  of  the  transition,  and  /  is  the 
oscillator  strength  of  the  transition.  The  rate  of  depopulation  of  level  n  due  to 
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photon  absorption  to  higher  energy  levels  is  equal  to 


dP_  N 


-  =  -  f  0.02654  Pn  0  /  sec-1  . 

„  u, ,  n  mn  nm 


m=n+l 


The  rate  of  population  of  level  n  due  to  photon  emission  by  higher  energy  levels 
is  equal  to 


dP  N  . 

E  P  Asec-1  . 
dt  m=n+l  m  mn 


In  Eq.  (3),  Amn  (Einstein  radiative  transition  probability)  is  equal  to 

A  n  "  T1-  "•  66698  X  1016  /  X-2  . 

mn  g  '  nm  mn 

°m 


The  rate  of  depopulation  of  level  n  due  to  photon  emission  to  lower  energy  levels 
is  equal  to 


dP  n-1 

=  -  T  P  Anm  . 
dt  j-i..  n  nm 

m  =  1 


At  equilibrium  the  arithmetic  sum  of  the  rates  of  population  and  depopulation  due 
to  the  above  lour  processes  is,  of  course,  zero  and  we  obtain  for  each  of  the 
N  ..evels  an  equation  of  the  form: 


n-i  ,  . 

Y"  0.  02654  P  0  /  -PA 

'  m  nm  'ran  n  nm  > 


T  ( 0.  02654  P  0  f  -  PA  )  =  C 
,  V  n  mn  •'nm  m  mn  / 

m=n+l 


By  numerically  evaluating  the  coefficients  of  P  and  Pn  in  Eq.  (6),  we  obtain  a 
set  of  N  simultaneous  equations  of  the  form 


E.  anmPm  =  0  * 


(7) 


N 

If  we  include  the  additional  constraint  that  Y  P  *  1,  the  solution  is  obta..ied 

m 

m=i 

by  the  inversion  of  a  square  matrix  of  dimension  N-l. 

In  the  case  of  barium  ion  where  there  are  only  5  transitions  for  a  5-level 

term  scheme,  many  of  the  elements  are  zero  and  it  is  not  necessary  to  perform 

the  matrix  inversion  in  order  to  solve  for  equilibrium  conditions.  However,  for 

neutral  barium  there  are  62  possible  energy  levels  (if  we  include  the  ionization 

limit),  and  the  inversion  of  a  61  x  61  matrix  presents  problems  even  with  a 

large  computer.  In  addition  sir.ce  there  are  86  permitted  transitions,  all  but 
o 

172  of  the  61  elements  of  the  matrix  will  be  zero,  and  the  program  would  spend 

a  considerable  amount  of  computer  time  manipulating  quantities  of  zero  magnitude. 

In  addition,  because  of  the  interest  in  the  dynamic  properties  of  the  problem 

(the  time  it  takes  to  grow  in  to  equilibrium  of  the  ion  radiations,  the  time  constant 

1  3 

for  significant  pumping  of  the  '  D  metastable  levels  of  neutral  barium,  the 
possible  sensitivity  of  relative  intensities  to  the  identity  of  the  ion  precursor), 
it  seemed  that  a  stepwise  solution  to  the  problem  would  be  more  informative.  In 
this  approach,  it  is  initially  assumed  that  all  species  are  in  the  ground  state.  The 
sun  is  allowed  to  shine  for  a  small  time  6t,  during  which  a  certain  population  of 
excited  states  is  built  up,  depending  on  the  solar  irradiance  and  the  values  of  the 
oscillator  strengths  of  the  permitted  transitions.  Then  the  sun  is  "turned  off"  and 
the  excited  states  are  permitted  to  radiate,  thereby  adding  to  the  population  of 
certain  lower  levels.  By  properly  choosing  the  value  of  the  time  interval,  the 
result  will  be  a  good  approximation  to  the  situation  existing  at  the  end  of  the 
period  fit.  The  process  is  then  allowed  to  repeat,  using  the  new  relative  popula¬ 
tions  of  each  of  the  energy  levels  as  a  starting  point.  It  was  intuitively  felt  that 
the  correct  choice  of  the  value  of  fit  was  properly  made  when  it  had  a  value  some¬ 
what  greater  than  the  lifetime  of  the  excited  state  (that  is,  10  6  sec  for  weak 
transitions),  and  somewhat  less  than  the  lifetime  of  metastable  levels.  Although 
the  radiative  lifetime  of  metastable  levels  may  be  measured  in  seconds  or  tens 
of  seconds,  in  practice  the  actual  lifetime  if,  more  likely  to  be  of  the  order  of  one 

second  due  to  the  absorption  of  solar  photons. 

-3  -2 

Initially  fit  was  chosen  as  10  sec,  and  it  was  later  increased  to  10  sec 

without  significantly  changing  the  results  (when  examined  with  time  resolution  of 
0. 1  sec  or  greater). 

In  order  to  check  on  the  validity  of  the  method  a  simple  model,  which  could 
also  be  solved  analytically,  was  constructed  and  evaluated  b>  both  techniques. 

The  solution  of  a  2-level  1-transition  scheme  is  trivial  and  would  not  permit  an 
evaluation  of  the  technique.  Therefore,  a  3-level  2-transition  scheme  with  one 
metastable  level  was  chosen  for  analysis. 
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2.  ANALYSIS  OF  SIMPLE  3-LEVEL  SYSTEM 

2. 1  Exact  Analytical  Solution 

Figure  1  shows  the  simple  term 
scheme  and  permitted  tx-ansitions.  The 
solar  irradiances  were  chosen  so  that  the 
rate  of  population  of  level  3  from  level  1 
would  be  lx  Fj  sec  1 ,  and  the  rate  of 
upward  pumping  from  level  2  to  level  3 
was  0,2XP2  see'*.  In  addition  we 
allowed 

A32  =  108  sec"1 

A31  =  2  X  10B  sec1 

The  equations  governing  radiative 
equilibrium  of  this  situation  are 

f-  ■-5Pl+bP3 

dP9 

"dT  =  'c?2  +  dP3 


Figure  1.  Simple  Term  Diagram 
for  Hypothetical  Case 


8(a) 


8(b) 


dP„ 

HT  "  aPj  +  CP2  -  (b+d)P3 


and 


1 


where 


a  =  1.0 
b  =  2v  10° 
c  =  0.  2 
=  lv  108 


8(c) 


d 


Tat  fc-ssctl  at  ExpJ-  CVi  fcas  E*t«ai  pfrefcsaatdl ^  ts^atsada  flllliL  Psr  Ciw 

psrficc^r  css*1  wdssectftt  cccsSii5t2S2r«’altt*C'-  Pg  »  l„  P*  *P,g  *0*  tS» 

soJeearacs 

P,  ■  'j~,|i-fffit)j|  -*■«£)•*-  0»t=««9b£J 
ar-#r 

p_* sea 

p3*  -^sii-ffMj+xfc&j  5&; 

g~-a~ 

arbere 

t  2 

e“  -  ^  -  ac  -&•  ad  ■*•  be 


a~~  &  ■*■  b  ■*•  c  *  dj 

,3  -  (ac  *  *&  be) 

and  wfc-sre 

f(i)  -  c~sl  lco£:Jl+-sidJ>?l( 

1C&I 

g{0  =  e”*24  {crs&^!-jiid3,i) 

SXs) 

hM.e^-S^) 

** 

ir-fci 

V.  e  may  also  rewrite  Lq.  (10)  as 

«•»  -■?*"“  [("jK'M'-J X'] 

11(a) 

eft)  «}'’*  j|i-J ) <■••'* 

1  lift) 

h(t)  -•  “  c~at  [<*  -  •■'“t  j 

lift) 

11(c) 


* 


Thus 


(r^oc..j 


ZzzitzTUrx  veSazsfssrx,  &„  c  assdi  dS  me  tint! 

«  *  ILS*  11? -fr  CbS 
2_£r  12? ■*■  £T_  1233 

*  2.4  x  v?  -ft-  0.2 

fe  =  G.4x  rS® 

O? Cs£&  *  3  X  il?  -ft  C-2 

zrA 

Pj=r  Q.rtSii  -ft  e'°'-i5Sr'‘  jp.7143  -  0-023  V  ICf3] 

^e-3xlA  [a.1777xi0'S]. 

„g 

Oiroraiy  for  tunes  greater  Cx&rs  10  sec,  the  third  term  is  negligible  and  we  may 
write  tts  s-elut ion  as 


Pj  =  0-235*  +  e°-'5&S:€t  1 0.7143)  -  (12) 

12  V^wumir  mliiiol  >!>liDr 

It  appeared  that  an  approximate  analytical  solution  for  this  model  might  be 
possible  on  the  basts  of  the  following  argument: 

Of  the  atoms  excited  to  level  3,  two-thirds  will  decay  to  level  1,  and  one- 
third  to  level  2.  Then  Eqs.  (3)  become,  with  the  appi  opriate  numerical  co¬ 
efficients, 

dP 

at'-O  '  P2)  i/ » +  (o.  2  r2 )  2/3 


13(a) 


7 


dP_ 

~  =  -  (0. 2  P2)  2/3  +  <1 X  Pj)  1/3  .  13(b) 

Solving  tiwse  and  putting  the  initial  conditions  as  Pj  =  1,  P?  =  0,  we  obtain 

Pj  =  0.2857  +  0.7143  exp(-0.4666t)  14(a) 

P2  =  0.7143  [  1  -  exp(-0.4666t)]  ,  14(b) 

Note  that  the  solution  Eq.  (14a)  is  identical  with  that  obtained  as  Eq.  (12).  Thus 
this  approximation  is  valid  for  all  but  the  first  microsecond  or  so,  and  it  is 

therefore  possible  for:  (1)  the  solution  to  the  dynamic  grow-in  to  equilibrium  of 

the  barium  ion  to  be  derived  as  the  solution  of  3  (rather  than  5)  coupled  linear 
differential  equations;  and  for  (2)  the  corresponding  neutral  barium  process  to  be 
obtained  as  the  solution  of  a  set  of  8  (rather  than  61)  equations. 

The  function  Eq.  (14a)  has  been  evaluated  by  hand  (using  a  Wang  program¬ 
mable  calculator)  and  gives  excellent  agreement  (at  0.  1  sec  intervals)  with  the 
computed  stepwise  solution  (to  be  described  below). 

Note  that  the  time  constant  for  pumping  of  the  metastable  level  is  2.  143  sec, 
whereas  the  simple  approach  used  earlier  (Best  et  al,  1970)  to  give  the  rate  of 
pumping  of  Bal  1D  would  yield  an  expected  time  constant  of  3  sec.  Consequently, 
it  may  be  inferred  that  the  rate  estimated  by  Best  et  al  (1970)  for  Bal  1D  pumping 
may  be  inaccurate  by  as  much  as  30  percent;  that  is,  it  could  be  somewhat  less 
than  the  1. 24  sec  derived. 

2.3  Computed  Stepwise  Solution 

A  computer  program  has  been  written  for  the  general  ease  of  energy  levels  in 
the  range  l^N^NMAX  and  associated  optical  transitions  1<LsLMAX.  In  this  we 
assume  initially  that  all  the  population  resides  in  the  ground  state  N  =  1.  As  noted  in 
the  introduction  above,  the  computation  interval  DELTC  was  allowed  to  be  greater 
than  the  lifetime  of  the  excited  states,  but  less  than  the  time  between  optical  excita¬ 
tions  from  the  ground  or  metastable  levels.  This  model  was  allowed  to  be  pumped 
by  the  solar  irradiation  with  DELTC  =  0.  01  sec.  Since  this  was  a  test  of  the  general 
program,  a  brief  qualitative  description  is  given  here. 

Input  data  consists  of: 

NMAX  :  the  maximum  number  of  energy  levels 
LMAX  :  the  maximum  number  of  optical  transitions 
DELTC  :  the  computation  time  interval 
NMIN  :  the  number  of  stable  (or  metastable)  levels 


For  each  of  the  energy  levels  we  read  in: 

N  :  the  sequential  number 
E{N)  :  the  energy  in  cm  1 
G{N)  :  the  statistical  weight 

For  each  of  the  optical  transitions  we  read  in: 

L  :  the  sequential  number 

WVL(L)  :  the  wavelength  in  Angstroms 
NU(L)  :  the  upper  energy  level  number 
NL(L)  :  the  lower  energy  level  number 
F(L)  the  transition  oscillator  strength 
SINC(L)  :  the  solar  continuum  irradiar.ee 
FLD(L)  :  the  Fraunhofer  line  depth 

The  program  initially  checks  the  above  data  for  accuracy  by  computing  the  wave¬ 
length  of  each  transition  for  the  energy  level  numbers  of  the  upper  and  lower 
states,  and  then  comparing  with  the  input  data;  correction  to  wavelength  in  air  is 
made  above  2460 R  .  This  wavelength  is  used  so  as  to  permit  comparison  with 
the  data  tabulated  by  Miles  ct  al  (1969). 

Initial  populations  are  all  set  equal  to  zero,  with  the  exception  of  P(l)  =  1.  0. 
The  rales  of  population  of  each  of  the  upper  levels  from  each  of  the  lower  levels  are 
computed  in  turn  for  DELTC,  and  after  all  LMAX  transitions  have  been  con¬ 
sidered  corresponding  changes  in  upper  and  lower  energy  level  populations  are 

-2  - 1 

made.  If  II  is  the  solar  continuum  irradiance  in  watts  m  nm  and  is  the 
wavelength  in  microns,  then 

0  =  1679.4  X8  11  (15) 

and  the  rate  of  pumping  of  the  upper  level  NU  from  the  lower  level  NL  is 

dP/pj)  =  026g4  P(NL)0  f  FLD  .  (16) 

Then  the  system  is  permitted  to  radiate,  and  the  branching  ratios  are  computed 
in  terms  of  the  computed  Kinsfcin  A  coefficients: 

AE=G<nM0>  66698  y  lo8fXv‘2  •  U?) 


Since  all  lifetimes  are  much  less  than  DEI/l’C,  we  perform  a  check  to  ensure  that 
all  levels  with  N  >  NMIN  have  emptied.  The  emptying  process  for  each  level  is  in 
accordance  with  computed  branching  ratios.  Printout  and  plotting  of  data  is  called 
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for  at  meaningful  intervals.  This  data  includes  populations  of  all  stable  (or 
metastable)  levels  and  the  brightnesses  of  emission  lines. 

In  the  case  of  this  synthetic  model,  values  of  the  input  parameters  were 
chosen  to  match  the  input  data  of  paragraph  2.  J: 


N 

E(  N) 

GW) 

1 

0 

4 

2 

5000 

4 

3 

20, 000 

O 

6 

and 


_L 

WVL(L) 

F(L) 

SINC(L) 

FLD(L) 

1 

4998.60 

0.  375 

0.  959 

0.500 

2 

6664.  83 

0.333 

0.455 

0.  100 

The  computer  tabulated  output  agreed  exactly  with  the  analytical  function  computed 
according  to  Eq.  (14a). 


3.  It  \I1IB\I  ION  SYSTKW 


The  program  providing  the  stepwise  solution  to  the  optically  thin  radiative 
transport  problem  was  next  applied  to  the  barium  ion  system  with  5  levels  and 
5  transitions  (Figure  2).  This  allowed  the  equilibrium  value  to  be  compared  with 
that  calculated  earlier  (Best,  19G8),  as  a  further  check  on  the  program. 

At  equilibrium  the  values  obtained  were: 


N 

Etcm'1) 

P(N) 

P(N)* 

1 

0 

0.  575 

0.  575 

2 

4873. 85 

0.  130 

0.  130 

3 

5674. 82 

0.294 

0.  294 

4 

20261. 6 

5 

21952.4 

L 

WVL(A°) 

BR.(N) 

BIUN). 

1 

4554.03 

0.780 

0.780 

2 

4934.  09 

0.299 

0.299 

3 

5854. 08 

0.032 

0.  032 

4 

6141.72 

0.  243 

0.243 

5 

6490.90 

0.  103 

C.  103 

:  Values  (per  total  ion)  given  in  Best  (1908). 

+  Dr.  J.  H.  M.  Fu  (EG&G)  has  kindly  pointed  out  a  typographical  error  in 
Table  II  of  Best  (1908);  1. 1355  should  read  1,355. 
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This  equilibrium  was  reached  with  a 
time  constant  of  t  *  0. 869  sec,  with  the 
ground  state  population  decaying  as 


Kt  =  0.575 +0.425  exp  (18) 

i.  BARIUM  NEUTRAL  SI  STEM 

•1.1  Equilibrium  lilhoul  Pbjloionization 

The  program  was  now  run  with  input 
data  obtained  as  follows: 

(1)  Oscillator  strengths  were  taken 
from  the  NBS  publication  (Miles  et  al, 

1969). 

(2)  Solar  irradiance  S1NC(L)  was 
taken  as  the  value  for  the  continuum 
given  in  the  NASA  Atlas  (Arvesen  et  al, 

1969),  with  a  correction  FLD(L)  being 
included  for  the  Fraunhofer  line  depth  (relative  to  the  continuum).  The  FLD  (L) 
was  taken  from  the  Utrecht  Atlas  (Minnaert  et  al,  1940). 

Figure  3  shows  a  partial  term  scheme  for  Bal,  with  the  transitions  cor¬ 
responding  to  the  brightest  observed  emission  features  being  shown. 

Table  1  shov/s  the  input  values  used.  THETA(L)  is  the  continuum  solar  ir- 

-9  -1 

radiance  SINC(L)  expressed  in  photons  cm  "  Hz  ,  and  AE(L)  is  the  Einstein  A 
coefficient  derived  from  the  values  of  oscillator  strength  and  wavelength.  This 
model  was  computed  with  no  photoionization  process  included.  The  idea  was  that 
such  a  model  would  be  modified  by: 

(1)  Ilemoving  suddenly  one  half  of  the  metastable  population  (for  each  of 
several  metastable  levels  in  turn)  after  equilibrium  had  been  reached,  and  ob¬ 
serving  the  time  taken  to  return  to  an  equilibrium  situation. 

(2)  Allowing  photionizing  transitions  to  occur,  from  each  of  the  rnetastablc 
levels  in  turn,  with  a  rate  such  that  a  time  constant  for  the  process  of  about  30  sec 
was  obtained.  Then  the  relative  intensities  of  the  lines  for  each  possible  process 
could  be  examined.  It  is  expected  that  if  the  recovery  times  found  in  step  (1)  above 
are  short  compared  with  the  photoionization  time,  then  there  would  be  no 
detectable  difference.  The  photoionization  mechanisms  were  expressed  for 
computation  in  terms  of  an  oscillator  strength.  This  concept  is  valid  for  an  auto- 
ionizing  transition,  but  for  continuum  photoionization  the  rate  requires  to  be 
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Figure  2.  Barium  Ion  Term  Diagram 
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^iruro  3.  Barium  Neutral  Term  Diagram 
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Only  those  transitions  for  which  BR(L)  >  0. 00]  photons  sec  1  atom  1  are 
tabulated,  and  their  temporal  behavior  is  characterized  as  follows: 

Type  1:  Shows  "exponential"  decay  similar  to  P(l) 

Type  2:  Shows  transient  peak,  similar  to  P{5) 

Type  3:  Shows  exponential  grow-in,  similar  to  P(2,3  or  4) 

Type  4:  Falls  initially,  then  grows  into  steady  value 
The  rate  of  pumping  due  to  any  transition  is,  of  course, 

-  PWU  /nl>nu  theta  FLD  sec"1  (20) 

where  'ire2/mc  =  0.  02654  and  the  values  of  THETA  and  FLD  are  given  in  Table  1. 

1.2  Readjustment  of  Lquililirium  After  Sudden  Removal  of  Some  Mutnsluliio  Atoms 

Sirce  equilibrium  was  reached  within  7  sec,  in  this  run  one  half  the  (N=5) 

2 

population  was  removed  at  10  sec,  and  one  halj  the  Dg  (N=3)  population  at  20  sec. 
The  run  was  terminated  at  30  sec.  Recovery  from  the  first  perturbation  was 

complete  within  one  sec  for  the  ground  state  and  metastable  level,  and  within 

3  c  1 

about  2  sec  for  the  D  levels.  In  the  second  case  recovery  for  the  >D  level  was 

3 

also  complete  within  one  sec,  whereas  the  ground  state  as  well  as  the  D  levels 
required  about  2  sec  for  recovery  of  equilibrium  populations.  Since  these  recovery 
(that  is,  population)  times  are  so  short  compared  with  the  photo  ionization  time 
constant,  it  may  be  expected  that  the  relative  line  intensities  will  not  be  very 
sensitive  to  the  choice  of  metastable  level  assumed  to  be  the  precursor  ti  photo¬ 
ionization. 

1.3  IMioioionizniiun  From  ,,^l)  Mnt.istulile  Levels 

Photo  ionization  was  included  in  the  computation  process  by  the  following 
technique: 

(1)  Add  an  extra  energy  level  corresponding  to  the  ionization  limit 

(2)  Increase  NMAX  by  1  accordingly 

(3)  Allow  a  transition  from  the  lower  level  NL  concerned  to  the  ionization  limit, 
with  the  oscillator  strength  being  chosen  to  yield  a  time  constant  of  30  sec.  This 
hypothetical  transition  can  later  be  equated  with  a  photoionization  continuum 
process,  or  with  an  autoionizing  transition.  This  will  also  require  a  unity  increase 
in  LMAX. 

(4)  Permit  all  downward  transitions  as  before,  up  to  and  including  (LMAX-1) 


9033.99 
9215.52 
959G.  55 


15 


P(l)t  =  0.512  +  0.  394  exp  (-  £)  (19) 

where  the  time  constant  r  =  1.  067  sec. 

At  equilibrium  the  relative  populations  are: 


E(N) 

P(N) 

1 

0.0 

0.512 

2 

9033. 99 

0.  086 

3 

9215.52 

0.  141 

4 

0595. 55 

0.  176 

5 

11395.40 

0.  085 

These  have  been  reached  to  within  0.  001  within  7  sec. 

The  relative  emission  intensities  BR(L)  in  units  of  photons  per  second  per 
total  barium  atom  are  given  in  Table  2.  As  might  be  expected,  the  relative 
brightnesses  show  the  same  temporal  behavior  as  corresponding  energy  levels. 
The  5535A°Iine  shows  an  intensity  decay  similar  to  that  of  the  population  of  the 
ground  state. 


Table  2.  Equilibrium  Emission  Line  Intensities  per  Barium 
Atom,  Without  Ionization 


L 

WVL(L) 

Angstroms 

BR(L) 

photon/atom 

sec 

Type 

L 

WVL(L) 

Angstroms 

BR{L) 

photon/atom 

sec 

Type 

12 

2596 

0.  002 

1 

61 

6019 

0.  377 

3 

14 

2702 

0.001 

1 

62 

6063 

0.800 

3 

16 

2785 

0.001 

1 

63 

6111 

1.313 

3 

17 

3072 

0.071 

2 

64 

6342 

0.405 

3 

18 

3501 

0.  130 

2 

65 

6451 

0.  175 

3 

19 

3889 

0.012 

3 

66 

6482 

0.636 

2 

20 

3910 

0.021 

3 

67 

6499 

1.855 

3 

21 

3936 

0.006 

n 

O 

68 

6527 

0.  986 

3 

22 

3938 

0.004 

3 

69 

6595 

0.  828 

3 

23 

3993 

0.  144 

3 

70 

6675 

0.401 

3 

24 

3996 

0.012 

3 

71 

6694 

0. 4  84 

3 

25 

4132 

0.015 

3 

72 

6868 

0.  102 

3 

29 

4283 

0.077 

2 

73 

7060 

2.667 

3 

44 

4592 

0.066 

3 

74 

7120 

0.532 

3 

46 

4605 

0.  004 

3 

76 

7280 

1.  947 

3 

48 

4628 

0.  016 

3 

78 

7418 

0.064 

3 

52 

4726 

0.  124 

2 

79 

7488 

0.  367 

3 

53 

5535 

8.  112 

1 

80 

7672 

0.832 

3 

56 

5806 

0.057 

2 

81 

7780 

0.359 

3 

57 

5826 

0.441 

2 

83 

7911 

0.  128 

4 

58 

5907 

0.047 

3 

84 

15000 

0.  347 

1 

59 

5972 

0.368 

3 

85 

27750 

0.  048 

4 

60 

5997 

0.  383 

3 

86 

29223 

0.  124 

3 
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Only  those  transitions  for  which  BR(L)  >  0. 001  photons  sec  1  atom  1  are 
tabulated,  and  their  temporal  behavior  is  characterized  as  follows: 

Type  1:  Shows  "exponential"  decay  similar  to  P(l) 

Type  2:  Shows  transient  peak,  similar  to  P(5) 

Type  3:  Shows  exponential  grow-in,  similar  to  P(2,  3  or  4) 

Type  4:  Falls  initially,  then  grows  into  steady  value 
The  rate  of  pumping  due  to  any  transition  is,  of  course, 

-  4^”  4^=  PINU  ^  /NL>NU  THETA  FLO  sec’1  (20) 

where  ire2/mc  =  0,  02654  and  the  values  of  THETA  and  FLD  are  given  in  Table  1. 

4.2  ItuHiljuslMunl  of  Equilibrium  After  Sudden  Ilcmovnl  of  Sonic  Metnstuliie  Atoms 

Sirce  equilibrium  was  reached  within  7  sec.  in  this  run  one  half  the  (N=5) 

3  “ 

population  was  removed  at  10  sec,  and  one  hal£  the  (N=3)  population  at  20  sec. 
The  run  was  terminated  at  30  sec.  Recovery  from  the  first  perturbation  was 

complete  within  one  sec  for  the  ground  state  and  lD9  metastable  level,  and  within 

3  “  i 

about  2  sec  for  the  D  levels.  In  the  second  case  recovery  for  the  >D  level  was 

3 

also  complete  within  one  sec,  whereas  the  ground  state  as  well  as  the  D  levels 
required  about  2  sec  for  recovery  of  equilibrium  populations.  Since  these  recovery 
(that  is,  population)  times  are  so  short  compared  with  the  photoionization  time 
constant,  it  may  be  expected  that  the  relative  line  intensities  will  not  be  very 
sensitive  to  the  choice  of  metastable  level  assumed  to  be  the  precursor  ti  photo¬ 
ionization. 

hit  IMioloiuni'/.ntion  I'rom  ^1)  Mct.islablc  Levels 

Photoionization  was  included  in  the  computation  process  by  the  following 
technique: 

(1)  Add  an  extra  energy  level  corresponding  to  the  ionization  limit 

(2)  Increase  NMAX  by  1  accordingly 

(3)  Allow  a  transition  from  the  lower  level  NL  concerned  to  the  ionization  limit, 
with  the  oscillator  strength  being  chosen  to  yield  a  time  constant  of  30  sec.  This 
hypothetical  transition  can  later  be  equated  with  a  photoionization  continuum 
process,  or  with  an  autoionizing  transition.  This  will  also  require  a  unity  increase 
in  LMAX. 

(4)  Permit  all  downward  transitions  as  before,  up  to  and  including  (LMAX-1) 


The  oscillator  strength  was  initially  chosen  so  that 


"  dp'ml)' =  02654  F  THETA  =  30  sec_1*  (21) 

However,  it  soon  became  apparent  that  F  would  require  to  be  increased  by  l/P(NL) 
in  order  to  make  the  overall  time  constant  of  the  order  of  30  sec. 


Table  3.  Data  Yielding  Metastable  Photoionization  Time  Constant  of  Required 
Magnitude 


Precursor 

3d0 

2 

X 

X 

N 

3 

4 

5 

f 

0.  345 

0.  221 

0.  197 

AE  (sec-1) 

1.24(9) 

1.09(9) 

0.  619(9) 

WVL  (R) 

3046.33 

3082. 12 

3263.09 

SINC  (watts  m  2nm  J) 

0.539 

0.650 

1.  149 

THETA  (phot  cm-2sec-1Hz-1) 

25.59 

31.961 

67.  044 

FLD 

1.000 

1.000 

1.000 

t  (sec) 

34 

33 

34.5 

At  20  sec: 

P  (1) 

0.295 

0.295 

0.  284 

P((2) 

0.  04  G 

0.  046 

0.  049 

P'(S) 

0.  074 

0.075 

0.  080 

P  (4) 

0.  093 

0.  092 

0.  100 

P  (5) 

0.  048 

0.  048 

0.  04  6 

P  (G2) 

0.444 

0.444 

0.441 

Table  3  nows  the  values  of  oscillator  strength  required  to  achieve  a  photo- 
ionization  time  constant  of  about  34  sec.  Note  that  at  20  sec  the  relative  values 
of  P(N)  arc  practically  indistinguishable.  The  rate  of  ionization  for  a  bound- 
bound  (or  autoionizing)  transition  is  given  by 


1 

P(NL) 


dP(NL) 
dt 


=  C.02G54  F  THETA  FLD  sec' 


(22) 


18 


or,  if  it  is  a  continuum  phc  toionizing  process  the  rate  would  be 

1  l#n 


1  dP(N) 
P?NT  dt 


■  / 


\=0 


a  O^dk  sec 


-1 


j: 


aO  dv  sec 
v 


(23) 


where  a  is  the  photoionization  cross  section,  and  and  v^,  are  the  threshold 
wavelengths  and  frequencies  respectively. 

The  values  of  Einstein  A  coefficient  (AE)  corresponding  to  the  required 
oscillator  strengths  are  not  unduly  large.  However,  the  failure  by  Carton  (1962) 
to  detect  any  autoionizing  transitions  from  the  metastable  levels  renders  this 
process  unlikely.  It  is  therefore  more  appropriate  to  consider  in  detail  the  con¬ 
clusions  in  terms  of  continuum -photoionization  cross  sections. 

We  can  obtain  an  idea  of  the  cross  section  involved  by  making  either  of  two 
assun  ptions: 

(1)  That  the  photoionization  cross  section  varies  slowly  compared  to  the 
solar  irradiance  (which  falls  by  a  factor  of  2  every  100  to  200  A). 

(2)  That  the  photoionization  cross  section  falls  off  much  more  rapidly  (with 
decreasing  wavelength)  than  the  solar  irradiance.  Then  only  the  value  at  the 
threshold  is  significant.  In  either  case  we  can  say 

0.  02G54  /  =  Jo  dv~aAv  (24) 

where  Ay  is  the  effective  spectral  width  of  the  process  determined  either  (a)  by 
the  solar  continuum  intensity  as  a  function  of  wavelength,  or  tb)  by  the  photo¬ 
ionization  cross  section  as  a  function  of  wavelength. 

If  we  take  AX.  =  200$  (=  6  x  1013  Hz  at  3127/1 ),  we  find  that 

a  =  4.  423  X  10"1C  f  cm2 
=  442  f  Mb 

where  one  megabarn  (Mb)  =  10  cr.i  . 

Therefore  if  we  permit  photoionization  to  proceed  from  any  one  of  the  four  possible 
metastablc  levels,  we  require  a  cross  section  of  about  00  Mb.  If  we  consider 
photoionization  from  all  four  metastable  levels,  we  require  a  cross  section  of 
about  20  Mb  and  this  is  more  in  keeping  with  values  for  photoionization  from  the 
ground  state  measured  by  Hudson  etal  (1970'.  However,  these  values  were  only 
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maintained  over  a  few  X  .  The  fact  that  photoionization  was  probably  from  all 
four  metastable  levels  was  also  derived  from  the  results  of  the  ultra-violet- 
screened  barium  release  previously  reported  by  Best  et  al  (1969)  and  Rot; .  nberg 
et  al  (1971). 
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Appendix  A 


The  following  printouts  describe  a  computer  program  that  calculates  photo¬ 
equilibrium  of  barium.  One  set  of  printouts  describes  Program  Relax  which 
calculates  the  excitation  and  emission  of  atoms  in  the  transient  and  steady  state; 
and  a  second  set  of  printouts  contains  Sample  Data  to  be  used. 
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PROGRAM  RELAX 


program  ^;.!i<tr.MPJT,ojr?jr,TAp;5=r'<>jr,rftp;s=ojT?uT>TAPE?9tTAt>£2< 
lTA?£3,ra5E9) 

THIS  PROGRAM  CALCULATES  THE  EXITATION  ANO  EMISSION  OF  ATOMS 
IN  HE  TRANSIENT  AND  STEADY  STATE. _ 

INTEGER  3 _ 

COMMON  EC  IOC  >,3(103)  ,3(1J0)  ,P=»  ( 133  > 

COMMON  THETA (lOUi ,AEC1D0> , 9R ( 100) , RAE (133 > , RAT (100) 

COMMON  WV_  C100)  ,NU(103)  ,NL(100)  ,- (103)  ,sIN0(lCIQ>  ,FLOUOO> 

_ COMMON _ TIM(230)  ,3RT(233)  ,  PT  ( 20  0 ) _ _  _ 

OIMENSION  f  -A3L ( 3 )  •  TITl  ( 3)  ,  XL  (2)  , YL ( 2)  A3EK3)  ,YLABEL(3) 

0IMEN3I0N  mX(lOl)  ,INOY(101)  ,I?L<12)  ,  nTL£(8),SCALE(10) ,ISCL(3) 

_ SET  CONSTANTS. 

IND=0  _ 

'<=  0 

<K  =  0 _ 

OELTC=0.31 

TIME=3.0 _ 

WRITE (3, 10) 

0 _ F ORMAT  (HI,//) 

_ NMIN  IS  LOWEST  VALUE  Or  NJ  WHICH  SH3UL3  BE  £M3TY  AFTER  EMISSION. 

NMIN=9 _ 

REA 0  IN  THE  NJ.M3ER  OF  ENERGY  LE VE-5,  T HE  NJN3ER  OF  POS SI3LE _ 

'  TRANSITIONS  ANO  THE  NJM8ER  OF  SECONDS  THE  REACTION  SHOULO  GO. 

RFAO( 5, 23)  NMAX,lMAX,FTIME 

0  FORMAT  (2I5.F3.0) _ _ 

OO  11  NsZ.NMAX 

_ P(N) =0.0 _ 

1  CONTINUE 

P(l)=1.0 _ 

READ  THE  ENERGY  lEVELS  ANO  THE  3  FACTOR. _ _ _ 

INPUT  E (N)  IN  CM- 1,  G  IS  STATISTICAL  WEIGHT. 

RcAO( 5,30)  (E(N) , 3  CM) , N=1,NMA7>  *  “  ~ 

0  FORMAT  (5X>~ 10.0 , 13) _ 

_ REAP  THE  WAVELENGTH  OF  THE  TRANSITIONS,  THE  J3PER  AND  LOWER  ENER31 

L E V E _ S }  j  a  I NC ,  AHO  '■L3. 

WVL  IS  AN33TROM5  IN  AIR  Ir  XL  VAC  ,3E.  24E0A. 

Wft.  IS  AN33  T  ROMS  IN  VAC  Ir  WLVA'C  .Li.  24  50  ft. 

_ 3INC  IS  S3.AR  INTENSITY  IN  CONTINJJH  IN  WATTS  H»»-2  NH-1 _ 

( APPL .  0  =  T.  9,  345  (1970)  ) 

_ FLO  IS  rRA  JNHQFE  R  LINE  0 E 3 T H  » _ 


FLDISUOfnPaa  ND  A3SOR3TION 


READ C5, 43)  CHVLC-_),N:)CL),Ml.a)  ,rCL),5rMCCL),F1_DCL), 
♦L=l,LtffX) 

FORMAT  (5X» -10.0,215, ;i3.0,2f 1C. Cl _ 


CALCULATE  THE  VAlJES  3=  THETA. 

THETA  IS  53L&R  Ir.RAOIANCE  A*  KEAS  EARtr*  DISTANCE  IN  PHOTONS  CM* 
SEC-1  HZ-1  IS  SOLAR  C3MTIMJJH. 


DO  21  L=1,lMAX 

I-HETA  CL)  =SIMCCL) *I679.4*C  CHVLCL)  •1.3E-G4)  **3) 
CONTINUE 


CHECK  THAT  HE  ENERSY  DIFFERENCES  ARE  CONSISTENT  HITH  THE 
HAVElENST.hS  HHICH  HERE  PREVIOUSLY  READ. 

DO  31  L=1,LMAX 
NUL=MU(L) 

NLL=NLCL) 

HLVAC=1.0E0 S /  CECH'JL)  -E l ML-_)  ) 

IFCHLVAC.Lr.2460.)  WLD=WL/AC-HVL( L) 

IFCHLVAC.SE. 2460.)  HL4R=C<iLVAC-f).  075) /i. 000265 
IF CHLVAC.3E. 2460.)  HL3=HLAR-HVLCL> 

IFCA3SCHlD) .ST.0.25)  S3  T3  5 

GO  TD  31  ~ 

HRITEC&t SO)  L,HVL(L),SUl,SLL,ECSJL),ECSLL) 

5  FORM  AT  l*  T-iE  FOLLOWING  DAIA  IS  NOT  CONSISTENT,*//, 

+*  L  =*I3,*  HVLCL)  =*F3.2,*  HUCL)  =*I3,*  NLC-)  =*I3, 

+*  ECNOCD)  =*F8Tl,»  ECML(-))  =*FS.l) 

IND=1 

EXIT  IF  ANY  VAlOTIS  nDT  CONSISTENT. 

TO1,  £D'.Tj“?T35 - 

Calc jlaI'e  havelensths  in  microns,  and  eINstein  COEFFICIENTS. 

“0jT91  L  =  1,LHAX 
NUL=NU(L) 

*  NLl=NL"C  l) 

WLM=1.0£*4/CECMJL)-E(NLL> ) 

'  AE(L)=CFLOATCGCNLL))/rLOAT  C3  CNUL)  ))*3.6&98E*8*=CL)/CHLH**2) 
CONTINUE 
HRITE (o, 110) 

FORMAT  C4X,*L*3X*HVLCU *3X*NU*2X»NL*6X*FCL) *3X 
+*SINC CL) *2X* THETA (L) *4X*F.D(t) *4X»AE(L) *)  “ 

HRITE <5 ,200)  (L ,HVLC L) , NJ (L) ,NL (■ ) , FC L) , 5INC (l) ,THETACL) , 
+FL0CL) , AECL) ,L=1,LMAX) 

_ FORMAT  CI5,;10.3,2I4,4F10.3,E10.3) 


** 

t » 

**■ 

:»1CJ.UE  HE  **£  23£?ir£2££t!g5  ME2-S  !>£«£*«;%£  *2*  T*E  t£«r£iS 
tans.  3E  3£»i£f£j  3*  g*X  £*15513*  F3M  H£  ILstfEt. 

«** 

ca  iai 

SCM*3.® 

33  ill 

in 

If(KJ(!:3.E1.1l  SJUSflHtAEliJ 
:3lffItJE 

3£C=5j1*3-.I; 

IF(3EC.lT.13.II.At3.K.3E.<TI1l}  W2IT2C3.UH  <S 

S3 

?»««  C/*  _£V£iJ£3*  MU  <831  SJM’J-fEU  £t»Il£0f» 

33  121 

121 

IFCKJILI.E3.tl  iEAECL>=*ECLS/5>J< 
caxiiwE 

131 

cctritJE 

3  3I*?J  l3D»  SHIH  x:j.*l£s  HE  2?  SfieES  I'D  HE  E3I5SI3V 

3  EYEHf  3E_fC  5£22t3S. 

jr 

55 

TI»E=TI1EOELTC 

03  at 

??CN) =3.3 

31  C3NTINJE 


<=<♦1 


CALCJL&TE  WE  Or  ?WIIG  AT  EA2S  L£/£L. 


-* 

33  41  _=1 » -MIX 

NJl=XUCl) 

NLL=tLCLI 

.3»Tai=0.32a5%*:>CNLt.i*st<-,»"»H£*«C.|,*-£.bC-»*3£i,<C 
??  CMJL) =>?CtJLJ  **ATCL» 

41 

P?tNLLl=^c*j_L>-2ATCU 

CONTINUE 

L'ALCJLATE  NE4  P3LJPATI3NS. 


_ aa  7i  n=i,nmax 

:»IM)=?(N)  *33CN> 

_ PP(MI =0.3 _ 

71  CONTINUE 


CALCULATE  EMISSIONS  F?3H  EACH  LEVEL. 


LHAX1=_NAX-1 
30  51  L=l> -MAXI 
NUL=NU(L) 

NL.=NL(.) _ 

3  3*=  =M0T0N5/SEC. 

SATE=PCNJL) *W(L) 


we.s.1 

PPtXUS  ~XKf 


>1  £*j= 


S3  51  X=1,««X 


I'lK.S-.XIi  t.  *.*J> 
33  13  51 


W?tI£io,3ll  * 

?S»*F  M*  .£JL.»I3»  tOI 


ir>vin 


IfliC.E3.il  33  33  55 
Ef«C/2®l”3.*E.i«l  30  fO  53 

_ tfm!C.-E.2.tl  33  FO  55 

<f(«/ieil*13I.S£.<l  33  13  55 

_ If<n!g..:.13.SI  33  T3  65 

55  33x1*531 

«=«*!  _ 


«!IE  3W  31  *  F*?E  =32  5F532A3E. 


KlfEiEl  THE,  132J11 


73 

— xziTEti,f£)  rise - 

FOZ11T  t//,58X,*F11E  =*?7-3/l 

133 

«iHTo^33i 

•3?«T  (/,  V^/VSX’EOO  *SX*Pf «  */» 

S3 

M3ll£<5»53Z  S1,E(«,?(11,1=i,F«1)»11«,£H1»X»  jfCKHlX) 

FOJfcHT  (53X,X5,2F1Q.3> 

1V3 

K2IT£t5,l<*3J 

F3211T  C//,3W,*L*^X*«Vi.ll»»8X*2»I£Ci.)*3X*XU*3X*WL»10Xi'a2tL)»/> 

K2IT£t5,l23J  tL,H*LCL»,S*r  Is.)  ,HLCL>  ,3*CLJ  ,L=l,LMAXj 

12e  F02HT  t33X,i5»Ft9.3,E15.3,2I5,r15.3> 


CH£C<  If  41  -UV£  aLL0££g  Trf£  CH3Jt.*TI0\  T3  S3  TO  COMPLETION. _ 

ir(Ti«;,.r.FriM£i  so  ra  55 _ _ 

EN3ITLE  2 

P.ZMIH3  2 _ 

NOTE:  The  remaining  portion  of  the  program  sorts  the  data  so  that  it  can 
be  plotted.  This  portion  of  the  program  is  particular  for  the  CDC6F00,  and 
subroutine  PIX)T  V  is  a  plotting  subroutine  written  by  the  author  for  the 
CalComp  Platter. 


set  jp  rass tints  "02  ?Lorms. 


^^f^>^EILff  W,'fl'*1  V^l  i»||^-i^^|  nimrTiyfpvwu1  *V'Y'^yft'MI,mKV’*  1  .tttt 


23 


3AIA  X_/3.3,30.3/ 

TLtI>=0.3 

XS=5.9 

3ATA  X.A3E./3QH  TINE  (SE3S) 

/ 

DATA  Y.A3E./30H  3RISHTMESS  / 

DATA  YLA3L/33H  P0PJLATI3N 

/ 

DATA  ITITLIIJ  ,1=2, 3, 2)  /4+10H  ISN-1)  / 

OATS  CTir*_£Cr>,I=2,3,2>/4*10HAN3STR3iS  / 

- 

DATA  EH/IOHWAVELE.NSTH  / 

•* 

OPEN  MASS  5T3RASE  FILE  ON  JNIT  8  r0R  ROTATIONS, 
FOR  3RI3HTME5SE3,  AND  STORE  THE  DATA  ON  TA  =>E  2  IN 

ANO 

THE 

ON  JNIT  9 
PROPER  FILES 

n* 

CAL.  3?£NiSt8,IM3X,131,0> 

CALL  0?ENN5t9,IN3Y,101,3> 

33  131  L=1,1NAX 

33  151  __=1, « 

REA0C2 1  THC.L)  ,(3INCi*t),r-0(M>  ,N=i,„MAX) 

3RULL)=5IN31L) 
?T (LL)  -rt-Oli.) 


151  33MTIMJE 

_ CALL  WRITiS(9,3RT, :«,u> 


IFTL.3T.NiAX)  SO  TO  135 

CA^L  HRirN3(S,PT,<K,LJ 

135 

REWIND  2 

121 

CONTINJE 

IN3=0 

DATA  SOA.E/2. 5,  1.0, 0.75,0. 50, 0.25, 0.125, 0.05,0.  025,0.010,0.005/ 

r* 

.« 

READ  HE  300E  FOR  THE  TY=>E  OF  DATA  TO  BE  PLOTTED, THE  NUMBER 

•# 

OF  S3E0TRA.  _INES  OR  ENERGY  LEVElS,  HE  300E  FOR  THE  SPECTRAL 

•s 

LINE  OR  ENERGY  LEVEL  ANO  HE  SCA^E  CODE. 

23 

REAOt  5,100)  ENG,  J,  (IP- (I)  ,1  =  1,4)  ,150.(1)  ,  IPlI  I  >,1  =  5,8)  , 

1IS0LC  2) ,  CPt(I)  ,1=9,12)  ,  ISOl  ( 3) 

100 

FORMAT  (A10, 110,1514) 

IS=0 

mJ 

«/ 

CHEC<  IF  A..  THE  OASES  HA/E  3EEN  =l0TTE0. 

** 

IF ( J. 01 . 12)  30  TO  75 

m* 

c 

CHEC<  I?  WE  ARE  GOING  TO  3L0T  ENERGY  LEVELS  OR  WAVELENGTHS. 

IU=8 

IF  (ENG. ED. EN)  Id=9 

00  141  1=1, J, 4 

N0=3 

14=1+3 

IF ( ISO. E3. 3 )  N0  =  i 

— 

-  —  --  . . 

L 


29 


IFdU.Ea.9)  90  TO  35 
IV=0 

DO  171  11=1,14 

TF^IFlTIT) - 

IV=IV*1 

IFdf.Ea.Q)  30  TO  35 
IF dl.GT .  J)  30  TO  85 


PREPARE  TITLE  AND  INITIALIZE  PLOTTER. 


IF(IJ.Ea.B)  ENCODEdO»l5a,TITL(2*IV-l))  E(I3) 

GO  T3  171 

85  TITL(2*IV-1)=10H 
171  CONTINUE 
IS=ISfl 
IS§=ISCl(I$) 

YL<2)=S.0*3C4L£d5S)-0.001 

rS=SCALE<I35l 

CALL  PLOTV{ND,XL,XS,YL,Y3,2,33,-l,-33,TITL  , XL ABEL, YLABL  ,8.5,11.) 
G0~nT45 
35  CONTINUE 
IV=0  ~ ’ 

00  151  11=1,14 
IP=lPL  ( 1  f) 

IV=IVU 

- IFTTTTjTrJJ  -j0-T5--55 - 

IF d3.  EQ.  0)  GO  TO  95 

- rFTTj;ET.^'_"Ni;isTEn-07t'5a7TrTimTrr-i:T)"'Hf-c-rP7 - 

150  FORMAT  (IX,FS.2,1X) 

SoTTTn 

95  TITLE(2*IV-1I =10R 

151  CONTINUE 
IS=IS+1 

- Tjs^rsirmm - 

YL(2)=8.Q*3GALEdSS>-0.001 
YS=Sd  Ale  (13  31 

CALL  ?L0TV(N3,XL,XS,Y-,Y3,2,33,-l,-3),TrrLE,X-ABEL,YLABEL,8.5,ll.) 


PLOT  CURVES. 


45  00  141  11=1,14 _ 

IFdl.GT.J)  GO  TO  141 

_ IFdPL(II).E3.0)  GO  T3  141 _ 

CALL  REAaH3(IJ,PT,KK,I?L(II) ) 

_ N0=2 _ 

N=II-I*1 

_ CALL  PLOTtf ( N3, T IM i XS  ,  PT,YS  , KK, N, 1 , -80 , TIT  L  , XL ABEL ,YL ABEL , 

18.5,11.) 

_ IN0=1 _ _ 

141  CONTINUE 

_ GOTO  25 _ _ 


75  CALL  PL0TV(4) 
STOP 

_ ENO _ 

HAR<  E'JCOUNTEREO 


52 

37  29. 

1 

2 

snrso  i 

9033.99  2 

'  3 

9215.52  a 

\  V 

3596*55  3 

4  $ 

6 

12266. C  1 

11  229 47.4 

12  2307V. 4 


pmi 


IV  24192.1 


15  2V531.5 

16  2V979.9 


17  25642.2 

18  25704.1 


19  25956.5 

20  26160.3 


21  26816.3 

22  28554.3 


23  30743.5 

24  30750.7 


25  30815.6 

26  30818.1 


9 


29  33905.3 

30  34370.8 


31  34493.9 

32  34602.8 


33-34616. 
34  34630.3 


3 


ItttI 


44.4 
38  35616.9 


39  35709. 

40  35762.2 


41  35785.3 

42  35892.5 

43  35894.3 

44  36200.4 


45  36495.8 

46  36628.9 


49  37774.5 

50  38499.9 


b  77 

7392.41 

20 

1  78 

7417.53 

12 

79 

7488.03 

11 

50 

7672.09 

10 

;  81 

7730.43 

10 

82 

7905.75 

20 

7911.33 

15000.00 


27750.22 

29222.72 


3263.09 


-01 

-01 


9.9*2-03 

1.002-01 


1.542-02 

2.632-02 

0.1570 


1.218 

1.178 


1.171 

0.297 


0.335 
0.330 
1.  149 


1.000 


.0 

.985 


00 


1.000 

1.900 


1.000 
J.,0  00 
l.'J'OO 


